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Abstract 

The orbifold generalization of the partition function, which would describe the gauge 
theory on the ALE space, is investigated from the combinatorial perspective. It is shown 
that the root of unity limit q exp(27ri/fc) of the q-deformed partition function plays a 
crucial role in the orbifold projection while the limit g — >■ 1 applies to M"'. Then starting 
from the combinatorial representation of the partition function, a new type of multi- 
matrix model is derived by considering its asymptotic behavior. It is also shown that 
Seiberg-Witten curve for the corresponding gauge theory arises from the spectral curve 
of this multi-matrix model. 
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1 Introduction 

The recent progress on the four dimensional M = 2 supersymmetric gauge theory reveals a 
remarkable relation to the two dimensional conformal field theory [1] . This relation provides 
the explicit interpretation for the partition function of the four dimensional gauge theory 
[21 [3] as the conformal block of the two dimensional Liouville field theory, and is naturally 
regarded as a consequence of the M-brane compactifications [H [5], which reproduces the 
results of the four dimensional gauge theory O [7]. It is originally considered in [1] for 
SU(2) theory, and extended to the higher rank gauge theory [8l[9], the non-conformal theory 
[TU1[TT1[T2] and the cases with the surface and loop operators [13], etc. 

According to this connection, established results on the two dimensional side can be 
reconsidered from the viewpoint of the four dimensional theory, and vice versa. One of 
the useful applications is the matrix model description of the supersymmetric gauge theory 
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[IllISllISllITllIHlliaEQiEl]. This is based on the fact that the conformal block on the sphere 
can be also regarded as the matrix integral, which is called the Dotsenko-Fateev integral 
representation |221I23) . In this direction some extensions of the matrix model description are 
performed by starting with the Liouville correlators on the higher genus Riemann surfaces 
[24:\ [25] . Furthermore another type of the matrix model is also investigated so far [26\ [271 128] 
[29]. This matrix model is directly derived from the combinatorial representation of the 
partition function by considering its asymptotic behavior while the Dotsenko-Fateev type 
matrix model is obtained from the Liouville correlator, so that the conformal symmetry is 
manifest. This treatment is quite analogous to the matrix integral representation of the 
combinatorial object Although they are apparently different from the Dotsenko-Fateev 
type, both types of the matrix model correctly reproduce the results of the four dimensional 
gauge theory, e.g. Seiberg-Witten curve. 

The purpose of this paper is to extend the remarkable connection between the two and 
four dimensional theory to the orbifold theory. The four dimensional orbifold manifold is 
given by C^/F where F is a finite subgroup of SU(2), e.g. F = Z^, and its minimal resolution 
of the singularity gives the ALE space j34j . The gauge theory on the ALE space is well 
investigated in |35[ [36] with respect to the instanton moduli space. Their results show the 
moduli spaces of the instanton on the ALE spaces are deeply connected to the affine Lie 
algebras, and the quiver varieties. This interesting fact provides a guideline to search for 
a counterpart of the AGT relation for the orbifold theory. Furthermore this theory has 
been reconsidered in terms of D-branes [37] in which some aspects of 2d/4d connection are 
partially studied, and the theory on the Taub-NUT space has been investigated in [38] in 
detail. 

To obtain the solutions of the four dimensional gauge theory on the ALE space, we have 
to deal with the partition function as the case of by implementing the ADHM construction 
and the localization method for the ALE space. By performing such a procedure, the orbifold 
generalization of the combinatorial partition function, which would describe the gauge theory 
on the ALE space, has been investigated [39] (see also [IQ]). They are defined as the invariant 
sector of the Young diagram under the orbifold action F, and the process to extract only 
the F-invariant sector is called the orbifold projection in [39J. This projection seems quite 
reasonable, but it is difficult to perform this in a systematic manner. In this paper we will 
show the orbifold projection is naturally performed when one appropriately parametrizes 
the partition function. This parametrization is just interpreted as the root of unity limit 
q — 7> exp(2Tri/k) of the g-deformed partition function while we have to take g — )■ 1 for 

^For example, the longest increasing subsequences in random permutations, the non-equilibrium stochastic 
model, so-called TASEP 30 , and so on (see also ;. Their remarkable connection to the Tracy- Widom 
distribution 32 can be understood from the viewpoint of the random matrix theory through the Robinson- 
Schensted-Knuth (RSK) correspondence (see e.g. [33] 1. 
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the usual theoryo The similar approach is found in [51] in the context of the spin 
Calogero-Sutherland model, where the combinatorial method plays an important role in 
characterization of the wavefunction, derivation of the dynamical correlation functions and 
so on (see, for example, [52]). 

In this paper we also propose a new matrix model description for the orbifold theory. 
Starting with the combinatorial representation and applying the method developed in |26t 
[271 EE] we derive a new kind of matrix models, and discuss the corresponding gauge theory 
consequences. As a result of orbifolding, we get A;-matrix models for the ALE space. This 
model is quite similar to the Chern-Simons matrix model [43j for the lens space S^/Z^ 
[g iSl He] (see also g^), which is recently applied to the ABJM theory [381 139] [SOl lST] . 
This is so reasonable because the ALE space goes to the lens space S'VZfe at infinity. We 
then discuss the large N limit of the matrix model, and show Seiberg-Witten curve is arising 
from the spectral curve of the matrix model. 

Along the recent interest in the 2d/4d relation, it is natural to search a two dimensional 
counterpart of the orbifold partition function. While the original relation relies on the 
interpretation of the SU(n) partition function as the conformal block of the conformal field 
theory described by Wn algebra, its generalization is also proposed [52l [53] [5l] [55]: the 
generalized W algebra corresponds to more complicated gauge theories in four dimensions. 
In such a theory an embedding p : SU(2) — ?> SU(n) plays an important role in characterizing 
the generalized W algebra [56]. On the other hand, in the orbifold theory, an embedding 
p : T C SU(2) — )• SU(n) is also found to characterize the decomposition into the irreducible 
representations of F. We will discuss this similarity from the viewpoint of string theory and 
propose a new kind of connection between 2d/4d theories. 

The organization of this paper is as follows. In section [2] we provide the combinatorial 
representation of the four dimensional gauge theory, which is originally proposed in [21 [3], 
as preliminaries for some extensions discussed in this paper. In section [3] we then study the 
orbifold generalization of the partition function with focusing on its combinatorial aspect. 
We will see the root of unity limit of the (^-deformed partition function implements the 
orbifold projection. Section U] is devoted to derivation of the matrix model from the orbifold 
partition function by considering its asymptotic behavior. In section [5] we study the multi- 
matrix model in detail by taking the large N limit. We then extract Seiberg-Witten curve 
of the corresponding gauge theory via the spectral curve of the matrix model. In section 
[6] we try to search a 2d/4d relation for the orbifold theory. The two dimensional theory is 
discussed from the viewpoint of string theory, and a relation to the generalized W algebra is 
also proposed. We also discuss topics related to the g-deformation of the partition function. 
We finally summarize our results in section [71 

^ Because the combinatorial partition function includes infinite products of the g-parameter, we have 
to take care of its convergence radius. Thus, to obtain the root of unity limit q — >■ exp(27ri/fc), we first 
parametrize it as g — >■ ujq, and take the limit of g — >■ 1. 
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Note added 

After submitting this article, some related papers, considering extensions of the AGT relation 
to the ALE spaces, appear in the preprint server [571 EH EH [60} l6T]. 



2 J\f = 2 partition functions 

Let us start with the partition functions for J\f = 2 theories, originally studied by [21 [3], 
as preliminaries of discussions along this paper: its deformation is proposed in section [3] 
and we will also see in section H) the matrix model can be obtained from the combinatorial 
representation by considering its asymptotic behavior. 

First we introduce the instanton part of the four dimensional partition function for M = 2 
SU(n) theory with Nf (anti)fundamental matters, 

^4d ^ ^A(--^/)l%i^,Z['l,,,,„ (2.1) 

A 

^4d ^ n r(Af ^ - ^ + + km + /3) r(/3(j -^) + h,) 

(M)iL) r(Af)-Ar)+/3(,-z) + M r(/5(,-.) + 6., + /3)' 

(2.2) 

antifund JU U U ^(6, - Mf + f3{i + 1)) ' ^ ' ' 

These combinatorial expressions are written in terms of n-tuple partitions, A = (A*-^-*, • • • , A^"^), 
and their parameters are related to those of the gauge theory as 

13 = - — , bi = — , aim = ai- am, km = k- bm, Mf = ^ (2.5) 
£2 £2 £2 

where ei > > £2 are fi-background parameters, ai is the Coulomb moduli parametrizing the 

U(l)"~^ vacua, and mj denotes the mass of the fundamental matter. We can also introduce 

bifundamental matter fields, but we do not focus on them in this paper. For J\f = 2* theory, 

which includes the adjoint matter, we have to consider another contribution, 

^4d ^ n r(Ay) - Af ^ + /3(j + km + m r(/3(j -^) + ku + m + P) 

(M)ik.) r(Af^ - Af) + /3(, -^)+km + Tl + P) r(/3(, -^) + k, + Tl) ^ ' ^ 



with the mass of the adjoint matter 
Note one can see a simple relation 



m = -. (2.7) 

£2 



^adj(m = 0) 
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This Af = 2* theory is given by the simple deformation of = 4 theory by applying the 
mass to the adjoint matter. We will comment on the massless limit of this J\f = 2* theory 
later. 

The five dimensional extensions of these partition functions are also proposed, which can 
be regarded as the g-deformed analogue of the original four dimensional one, 

^ "Z^^^ ^vec^(anti)fund' i^'^J 



1=1 f=ii=i [QiQnifq^ t »;g)oo 



-^fntifund — n ri ri _-, ' (2-12) 

1=1 f=ii=i {QiQmfq^^ t-'+^;q)oo 



(«,i)^(m,j) {QmQlmq * ' ^joo 

Here we define 

00 

{x;q)oo = ll{l-xqn, (2.14) 

p=0 

and deformed parameters are given by 

g = e^^ t = e-^i=/, Q; = e'^'=(7'S Q«™ = e"'"'^™ = (z'"", (2.15) 

Q^^=e"'f =q^f, Q, = e'" = 9^. (2.16) 

The radius R of the compactified dimension is implicitly included in these parameters 
by rescaling the parameters, e.g. q = e^"^^. Thus one can check the four dimensional result 
is reproduced by the five dimensional partition function by taking the limit ii —> 0, or 
equivalently g — t- 1. 

To obtain the matrix models, it is useful to consider another representation for these par- 
tition functions with cut-off parameters N^^^ for the partitions \26 \ 1271 [28]. By decomposing 
the partition functions ()2.10p . (I2.13P as 



n^vir\ (2.17) 



l,m 
n 



u = n^Sj ' (2-18) 

l,m 



5 



each part can be rewritten as 

^{1,1) 



"adj 



1 n (, 



t'-';s)=„ '?r' (<,^;"t'""-'+V 



n 



(2.19) 



n 



^(()_^(m) 



ii ^(0_x(0 ... ii 



Ar("i) 



jy(0 ^(m) 



,(«,m) 



nn 



{QxaQlmq 



1 (Qmg^''t^'"-^+^;g)oo 



(2.21) 



^(0_^(m) 



(2.22) 



We can check finite constants turn out to be 



A, 



{l,m) 



{Qmt; I 



iV(0 



(Qm;g)o 
1 

n iQrnQlmf;q)^' (iV«<ivM) 



(Z = m) 

(iVW =ivM,Z^m) 



n (<5mQimi-*+^g)oo (iV«>iVM) 



i=l 



(2.23) 



Therefore they can be represented as 

1 \ -i-r iQlmQ 



^(i)_^(m) 



n 



;,m ^m=0 / {l,i)j^{m,j) [Qlmq ' ' \ q)oo l,m i=l [Qmiq ' * 



n 



P-';q) 



n NW 



.(I) .(m) 



nn 



=0 TTTT (Qi™?*'" 9)» 



, (2.24) 



n^m-'M n 



x(0 x("i) . . , , , 



00 TT TT IVrnVmiy > yjoo 



(i,i)7^(mj) {QmQlmq 



. (0_,(m) 

i i 



nn 



g)oo i,m i=l {QmQlmq * ^ 



Af ^7V(™)-j+l. , 



(2.25) 



^^r. = nnn , '^'^X.rr . ■ 

/=i/=ij=i [QiQmfq^ ^ t ^■,q)oo 

^antifund = 1111 II ^- ,il ^^'^^^ 

i=i/=ii=i [QiQmjq i t '+^;g)oo 
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These expressions are convenient to give the matrix model description [26 \ \27\ [28], We will 
investigate deformed versions of the partition functions, starting from these expressions, and 
derive the corresponding matrix models in section HI 

3 Orbifold partition function 

We then consider the orbifold generalization of the partition function for TV = 2 theory 
[39] , which describes the gauge theory on the ALE spaces [Ml ESJ [36], obtained from the 
minimal resolution of orbifolds C^/F where F is a finite subgroup of SU(2). This resolution 
of the singularity is performed by replacing two-spheres, where their intersecting numbers 
are related to a Cartan matrix of the corresponding Lie algebras q. This connection between 
subgroups r of SU(2) and Lie algebras g is known as the McKay correspondence. Especially 
the Abelian groups L = correspond to A^^i Lie algebras, and we focus on them in this 
paper. 

To study such an extended partition function, we will show it is useful to deal with the 
(^-deformed partition function which corresponds to the five dimensional theory, and take 
the root of unity limit of the deformation parameter. We remark the procedure used in this 
paper to obtain the orbifold partition function is closely related to the method proposed 
in [H] to study the spin Calogero-Sutherland model (e.g. see i42j), where combinatorial 
arguments are quite important as well as the four or five dimensional gauge theory. 

3.1 Orbifold projection and root of unity limit 

To construct the orbifold partition function, let us first clarify the orbifold action for the 
combinatorial partition function. Recalling the partition function is obtained from the Chern 
character of the fixed point in the moduli space under the torus action, or the equivariant 
if-theory class, it is reasonable to observe an effect of the orbifolding on the torus action. 
Their expressions are closely related to the five dimensional partition function, so that we 
first study the five dimensional function rather than the four dimensional version to discuss 
the meaning of the orbifold action. 

By considering ADHM construction and the localization method for the ALE space 
[(i'2\ [5U] as in the case of the usual space, it is shown the orbifold partition function 
should be defined as a sector of the g-deformed function (with q ^ 1), which is invariant 
under the following orbifold action 

r : q — > ujq, t — > ojt, Qi — > uj^^Qu (3.1) 

where uj = exp{2TTi/k) is the primitive k^^ root of unity, pi is an integer pi G {0, • • • , A; — 1}, 
and pim = Pi — Pm- It denotes the 1*^* Chern class of the instanton bundle [39], and also is 
interpreted as the linking number of the five branes j38j . Actually the boundary of the ALE 
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Figure 1: T-invariant sector for U(l) theory with A = (8, 5, 5, 4, 2, 2, 2, 1). Numbers in boxes 
stand for their hook lengths h{i,j) = Xi — j + Xj — i + 1. Shaded boxes are invariant under 
the action of F = Z3. 

space is a lens space S^/T, and thus we can assign non-trivial flat connection at infinity in 
this case. Therefore the instanton solution is classified by both the 1*^* and 2"^^^ Chern classes. 
We remark the instanton number defined as k/\T\ coincides with the 2^^^^ Chern class only if 
the F^* Chern class is vanishing. This constraint is quite non-trivial, but it does not concern 
our derivation of the matrix model discussed later. 

It is easy to see this action in terms of a Young diagram. Since the five dimensional 
partition function is written in the following form 

n. 00 ^ 
l,m i,j 1 - QlmQ^ :> 

where Xj stands for the transposed partition, we can see the F-invariant sector satisfies 

Xf^ + X'f^ -i-j + l+pi,r, = (mod A:). (3.3) 

For U(l) case the left hand side coincides with the hook length of the box, h{i,j) = Xi— j + 
Xj — i + 1. Thus the corresponding four dimensional function can be given by taking into 
account this condition, 

i,m r-inv. aim + ei[t - X) ' -1) + e2{Xl - j) 

The product in this expression is explicitly taken over the F-invariant sector ()3.3p . Fig. [1] 
shows F-invariant sector for U(l) theory with F = Z3. 

Although this definition seems natural for the orbifold version of the partition function, 
it is not useful anyway because we have to extract the F-invariant sector by hands. This 
procedure is called the orbifold projection in [39]. On the other hand, when we consider 
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another deformation of the partition function, which is given by replacing the parameters 
with (|3.ip as 

n oo ^ 

■^vec ~ II I I -{0 , !(m) '. . , , , ~m ?(m) ~~r' (3-5) 

«,m i,j 1 - CJ ' ^ J QlmQ ' ■'t J 

and take the hmit of g — ?> 1, namely the root of unity limit q — )■ exp(27ri/A:) of the original 
(/-partition function, we can see only the F-invariant sector contributes to the partition 
function and the others are decoupled in this limit. It is because, unless the power of uj 
becomes xf^ + X^J^^ —i — j + ^ + Pim = (mod k) in (j3.5p . they just give a factor independent 
of the shape of the Young diagram, or the partition. Thus, if we take into account the 
adjoint matter contribution to regularize the singular behavior at g ^ 1, the weight function 
behaves as 

^ a;^+ei(»-A^"^'-l)+e2(Af>-j)+m (Z) - (m) 



J 



1 if \f + A^"^ - i - j + 1 + p/™ ^ (mod A:) 

(3.6) 

This means the orbifold projection is automatically assigned by this parametrization. There- 
fore let us define the partition function modified with ()3.ip as the ((/-deformed) orbifold par- 
tition function from now. We note that the pure Yang-Mills contribution can be extracted 
by taking the decoupling limit m — )■ oo. 

We now check this reduction with a simple example, SU(2) gauge theory on C^/Z2 with 
the adjoint matter. If we set €2 = — ei = h, ai = —02 = a and pi = P2, lower degree 
parts of the instanton partition function are obtained by replacing (/-parameters as q ^ —q, 
Q12 Q12 [39], 

^ (1 + Qmg)(l + Qmg-^) (1 - QmQl2)(l - QmQ2l) . 
° (l + g)(l + g-l) (1-Ql2)(l-Q2l) ' ^'^ 



(1 + QmQl2q){l + QmQ2iq'^) (1 + Qm<3l29"^)(l + QmQ2iq) / (1 + Qmg)(l + QmQ'^ 
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~ {I + Ql2q){l + Q2iq-^) {1 + Ql2q'^){l + Q2iq) V (1 + (/)(! + (/"l) 

^ ^ (1 + Qmg)(i + Q^q-') (1 - Qmg')(i - Qmq-^) 

™ (l + 5)(l+g-l) (l_g2)(i_g-2) 

^, (1 - QmQl2)(l - QmQ2l) (1 + QmQl2g)(l + QmQ2ir^) 
(1-Ql2)(l-Q2l) {1 + Ql2q){l + Q2iq-^) ' 

Thus, taking the limit g — )• 1, we have 

m2 



(3.8) 



(3.9) 



^ (3-10) 

2 \ / 2 \ 
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a; = exp ^ 



\1 




1 47r 

cj^ = exp =1 



Figure 2: Parametrization of q corresponding to = 2 theories on and 'C? jTL-^ with the 
Umit of |g| — )• 1. |g| is related to the radius of the compactified dimension . 

These correctly reproduce the result derived in [39]. 

Let us consider the meaning of this parametrization without taking the four dimensional 
limit \q\ 1. A candidate of the corresponding manifold is the Taub-NUT manifold because 
it can be obtained by compactifying the singularity on S^, and reproduces the ALE space by 
taking a certain limit. It is just a speculation, thus we have to come back to this identification 
problem in a future work. 

We then comment on ambiguity of this parametrization. In this paper we apply the 
primitive root of unity to define the partition function as ()3.5p . but other roots of unity, 
represented as w'" = exp(2r7ri//c), are also valid for the orbifold projection while k and r are 
co-prime. This arbitrariness reflects symmetry of the system. 

Fig. [2] shows how to parametrize q to obtain the partition functions for the ALE spaces. 
\q\ corresponds to the radius of the compactified dimension S^. When we consider the four 
dimensional limit, we have to approach \q\ 1 with taking care of the radius of convergence 
since the partition function includes infinite products of q. In [39j a similar parametrization 
is actually proposed, which is interpreted as an analytic continuated one, but they do not 
take into account regularizing the infinite product appearing in the partition function. Thus 
the parametrization used in this paper should be more suitable. 

3.2 Rearranging partitions 

As discussed before the partition function for SU(n) theory on is represented with n- 
tuple partition. For this case, we can reproduce the one-matrix model by blending n-tuple 
partition to a single one \26\ \T7\ I28j . Next we consider its natural generalization to the 
orbifold theory. In this case, on the other hand, it is convenient for the discussion below to 
divide a n-tuple partition into a /cn-tuple one [30], 
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Figure 3: The decomposition of the partition for A; = 3. 



where 



fc-i 

r=0 



(3.13) 
(3.14) 



r=0 



This corresponds to the decomposition into irreducible representations of F. 

Let us practice with an example for U(l) theory as shown in Fig. [3l Starting with a 
partition A = (8, 5, 5, 4, 3, 2, 2, 2, 1), we obtain a configuration (15, 11, 10, 8, 5, 4, 3, 1) via the 
mapping Aj \i + N — i + p with = 8 and p = 0. This is interpreted as the mapping of a 
Young diagram to a Maya diagram, or equivalently bosonic to fermionic variables. Classifying 
entries by modulo 3, we have three configurations (15,3) — > (5,1), (10,4,1) — )■ (3,1,0) and 
(11, 8, 5) (3, 2, 1). We then reproduce three Young diagrams from them by specifying A^^*") 
and p'^'^^ . 

We define the explicit relation between elements from both sides of (j3.13p as 

i + Ara,-)+/,0)+^^A?)-i + iVW+pz with j = cf'\ 



(3.15) 



Here c\ ' means the mapping from the index of the divided fcn-partition to that of the 
original n-partition. Introducing another set of variables 



i + 



Ar(/,0+/,r))^ /if ^aJ 



.(0 



i + 



+ Ph 



bi-pi+ r, 



(3.16) 
(3.17) 



and setting N^^^ = N and iV^'-'^) = iV for ah / = 1, • • • , n and (/, r) = (1, 0), • • • , (n, /c - 1) for 
simplicity, the partition function (|2.17p is, up to constants, rewritten as. 



n / /i'''-/i'-'' + l ?i'''-/i<^''+(/3-l)(i-i)+6,,„-p,„+/3 ~x 
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The former expression is in terms of the original n-tuple representation, and the latter is 
written with a /cn-tuple partition. Blending the /cn-tuple to A;-tuple partitions as 



p(-r) _ ( nin,r) p{n,r) Jl,r) Jl,r) 

we finally obtain an expression in terms of A;-tuple partition, 

_ T-r [UJ' "q^ ^ ^^'^ J ' '■,q)c 

^vec — J_ J_ 



(3.19) 



Again c- stands for the mapping from the index of the /c-tuple partition to that of n-tuple 
one. We will discuss its matrix model description from the expression of (|3.20p in section HI 
We can also get explicit representations for matter parts in a similar way, 



00 



(^r-p,+1^4'''-^>.+»t+(/3-l)(N-cW)+/3. 



i=l r=0 i= 

4 Matrix model description 

We then derive matrix models from the orbifold partition functions. According to the result 
that the orbifold partition function can be written with a /c-tuple partition, one can see the 
fc-matrix model is naturally arising from the combinatorial expression. The multi-matrix 
model, derived in this section, is 

1 v-fe-iv-iv 



VXe =2 ^"-=0 ^«=i '■^^ (4.1) 
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fc-1 N , (r) 
r=Oi=l 

= Vvec(a;) + V(anti)fund(2;)- (4-3) 

We remark this matrix model corresponds to only the contribution from the instanton part, 
thus we have to introduce the perturbative piece when we consider the whole prepotential 
of the gauge theory. The deformed version of the Vandermonde determinant appearing in 
the matrix measure is given by 






with 



2 Sinh ^ ' ' ] nn 2sinh ^ ' ^ ' ] (4.4) 

/ r<s ij \ 

= ^(r-s). (4.5) 
k 

Here we introduce another parameter 7, which is related to /3 as /3 = A;7 + 1. Note that N 
stands for the matrix size here. This matrix measure coincides with that of the Chern-Simons 
matrix model on the lens space S^/Z^ |441 B5l H6] when we take 7 = 0, corresponding to 
/3 = 1. In section [5] we will study this case in detail. 

To obtain a matrix model from the partition function, we consider asymptotic behavior 
of the combinatorial expressions by introducing the following variables 

) = (4.6) 

£2 

and taking the limit of £2 — > 0. We study the asymptotics of the orbifold partition functions 
in the following. 

4.1 Matrix measure 

Let us start with the measure part of the multi-matrix model, coming from the combinatorial 
expression of ()3.20p . In this derivation we assume 

/3 = A;7 + 1 = 1 (mod /c), (4.7) 

in order to satisfy the condition 

= {ujqf. (4.8) 

This parametrization is proposed in thus we call it the Uglov conditior^ One can see it 
is quite important to obtain the matrix model description from the combinatorial expression 
because if we denote q = ujq, this condition implies 

{q,t)^iuq,ut) = {q,q^). (4.9) 



^ The reason why we assign this condition (|4.7p here is just a technical one. We will see, in the forthcoming 
paper ^63^, that this condition is not essential for deriving the matrix model description. 
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Due to the condition ()4.8p . the infinite product in ()3.20p is reduced as 
J-j- Q 



-^-^ / r-s+1 ^<'''-£<"^+(/3-l)(c<"'-c(''')+/3 ~^ 

= n n(l-a;'--^+V'^'-4^'+(^-^)('=^'^'-'''')+^) (4.10) 

Taking the limit €2^0 and using (|4.7p . we can treat them as 

1 - u^-^A'-^f'^ n (l - uf-^+^A'-'f''^ = (l - uf-^A-^'f'^j ( 1 - e'^^'-T'^ 

(4.11) 

where we use the identity 

fe— 1 

ll{l-u;^z) = {l-z''). (4.12) 

r=0 

As a result we obtain the deformed Vandermonde determinant A^(x) presented in (j4.4p . 

When we take the four dimensional limit, a factor including d^'^''^^ does not contribute in 
the leading order because it is expanded as 



2sinh I | = 2i sin ( — - — vr 



Thus the matrix measure goes to 



l-^l'^^^^Uotfll^vr 



(4.13) 



k—1 nN k—1 nN 

^'(-) ^ n n {-^' - n n {-^' - ■ ^'-'^^ 

r=0 i<j r<s i,j 

This shows there is no interaction between k matrices in this part when we take (3 = 1 
(7 = 0), corresponding to the most desirable SU(n) theory. We remark that it can be 
apparently written in a simple form in terms of /c-matrix model, but it is still difficult to 
write it down with the original combinatorial representation. 

4.2 Matrix potentials 

We then discuss the potential terms for the matrix model. To obtain the matrix potentials, we 
have to evaluate the asymptotic behavior of the quantum dilogarithm function [6l 



While the series expansion with q = e'''^ has been already investigated, an explicit expression 
for expansion around the root of unity, which is required for studying the orbifold partition 
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function, is not revealed yet. Substituting q = cje^^ into ()4.15p . we obtain a similar expression 

k—l oo 



log g{z;q) 



1 

€2 



r=0 n=0 

^Li2 ) + 0(62) 



m=0 



(A; - r)"-"^F 
(n — m)!m! 



(4.16) 



Here Bn stands for the n Bernoulli number and Li„(z) is the poly logarithm function 



(4.17) 



For n = 2 it is especially called the dilogarithm. Since the expression of ()4.17p is valid where 
\z\ < 1, it is convenient to introduce the useful identity, called the inversion formula, 

^2 



1 vr' 
U2{z) +Li2(l/2;) = --(logz)^ ^ Y ~ 



(4.18) 



Thus we can extend its domain to the whole complex plane. 

Utilizing the results obtained above, we can evaluate the factor contributing to the matrix 
potential. The vector multiplet part in (j3.20p yields 



nnn 



(^r-p,+l^£W-6,+(/3-l)(N-c('-')+/3. 



- exp5:5:-iKt^,(xr)) 



r=0 1=1 



£2 



(4.19) 



Kec(^) = E [Li2(e'=(^-'^')) - Li2(e-'=(--"'))] + ©(e^) 

1=1 



n 



+ 1:2 5^Li2(e 



-k{x-ai) 



)• 



(4.20) 



Here we neglect a redundant constant term in this potential, and if x — < 0, we have to 
redefine this by using the identity (|4.18p . The derivative of this potential is given by 



1=1 



2sinh ( -{x - ai) 



(4.21) 



We can also derive the potential for the four dimensional theory. Taking the four dimensional 
limit, we have 



2 



and its derivative 



Kfc'{x) = lj2^og{x-ai). 



(4.22) 



(4.23) 



1=1 



We will show this potential plays an important role in obtaining Seiberg-Witten curve as the 
spectral curve of the matrix model. 
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We then investigate the matrix model potentials for the fundamental matter from the 
expression ()3.22p in a similar way, 



k-l Nf nN , i+,.-{Mf+l-f)i).^\ k-lnNo 



r,(, ,-l+r„-(r' '+Mf + (l3-l)i+l). z\ - i-^ i-^ g2 



(4.24) 



^fund = ^X;Li2(e'=(-+'"/)), (4.25) 



N 



1 ^ 

^fund = -- J^[(x + m/)log(x + m/)-(x + m/)]. (4.26) 
/=i 

We again neglect a finite constant independent of x for them. For k = 1 these potential 
functions coincide with the usual matrix model potentials for Af = 2 theories |261 128) . We 
now remark this matrix potential is independent of the deformation parameter /3, and its 
dependence only appears in the matrix measure. 

4.3 TV = 2* theory 

Let us derive another matrix model by taking into account the contribution from the adjoint 
matter, which should describe A/" = 2* theory [27] . After almost the same procedure discussed 
above, we now obtain the following matrix model: 

Z = J VXe~^^'^'-='>^^='^'''''''^\ (4.27) 

r=Oi=l ""^ ' 

^(x) = Kcc(^) + Kdj(x). (4.29) 

Again stands for the matrix size here. In this case the matrix measure requires another 
contribution, 

Ai(x) = nn 2sinh^ ^ nn 2sinh^ 

t:} ^ ( k{xt^ - xf ) + m) V Jl / fe(xf ) - xf + m) V 

r=0 i^j \ J r^s t,j \ / 

(4.30) 

The matrix potential for the adjoint matter is 

„ n 

K.dj(^) = -^(x + m)2 -^Y.^h (e-'=(--«'+-)) . (4.31) 

1=1 
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As a result, the total potential function becomes 

2 

1=1 



„ n 

Kec(x) + Kdj(x) = -nmx + ^Y1 [Li2 (e^^'^^"'^')) - U2 (e-'=(^-'^'+-))] . (4.32) 



This result is consistent with [27] for k = 1. 

We now comment on a relation to A/" = 4 theory on the ALE spaces [40J. Since J\f = 4 
theory is obtained from J\f = 2* theory by taking the massless limit of the adjoint matter 
m — >• 0, the result of [40j should be related to the result discussed above. In the massless 
limit the matrix measure and the potential become trivial 



A2(x) 



1, V{x) — > 0. (4.33) 



This means, when we go back to the combinatorial representation, the combinatorial weight, 
represented in terms of the hook length, becomes trivial in the partition function. As a result 
it remains only the counting parameter, which is related to the dynamical scale. 

For U(l) theory, therefore it simply corresponds to the following partition function only 
with implementing the orbifold projection, 

Z ~ ^#{r-invariant boxes} _ (4.34) 
A 

This is similar to the partition function discussed in HQ], but not identical. We should 
investigate this type of partition function and clarify the relation to M = A theory in a 
future work. 



5 Spectral curve of the matrix model 

To study a connection between the gauge theory and the matrix model, in this section we 
consider the multi-matrix model in detail, which is defined as 



r=0 1=1 

where A'^{x) is defined in (|4.4|) . In this section €2 is replaced with Qs, and we focus on the 
case of /3 = 1. For a while we consider a generic potential V{x). The method used in this 
section is partially based on that developed for the lens space Chern-Simons matrix model 

5.1 Large limit and saddle point equation 

We are interested in the 't Hooft limit of this matrix model, in which 

gs^O, Af ^ 00, (5.2) 
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with fixing the 't Hooft couphng 

T = gsN. (5.3) 

Actually, in this large limit, the evaluation of the matrix integral reduces to the calculation 
of the critical points. 

If we define the prepotential for the matrix model as 

fc-l N 



r=0 i=l 

fc-l N ( ^(r) _ ^(r)- 



+2 log sinh 

r=0 i<i \ / 

+2 E E ^i^^ - — \ ' (5-4) 

which corresponds to the genus zero part, we would extract the information about the gauge 
theory, e.g. Seiberg-Witten curve, from this function. We can obtain the condition for 
criticality by differentiating the prepotential. 



= E ' ^ ' + E E ^ . (5.5) 

This saddle point equation is also given by the extremal condition for the effective potential 
defined as 

\ / s(7^r)i=l V / 

(5.6) 

This potential involves a logarithmic Coulomb repulsion between eigenvalues. If the 't Hooft 
coupling is small, the potential term dominates the Coulomb interaction and eigenvalues 
concentrate on extrema of the potential V'{x) = 0. On the other hand, as the coupling gets 
bigger, the eigenvalue distribution is extended. 

To deal with this situation, we now define the densities of eigenvalues for each matrix, 

1 ^ 

P^'\x) = -Y.^{x-xt^) (5.7) 

1=1 

where xi is the solution of the criticality condition (j5.5p . In the large N limit, it is natural 
to think these eigenvalue distributions are smeared, and become continuous functions. Fur- 
thermore, as in the case of the usual matrix model, we assume the eigenvalues are distributed 
around the critical points of the potential V{x) as linear segments. In particular there are 
n critical points for SU(n) theory. Thus we generically denote the l^^ segment for p^^^ (x) as 
C^^'^\ and the total number of eigenvalues N splits into n integers for these segments, 



iV = ^iV(^'^) (5.8) 



1=1 
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where N^'"''^^ is the number of eigenvalues in the interval C^^'^\ The density of eigenvalues 
p^^\x) takes non-zero value only on the segment C^^'^\ and is normalized as 

Ar('.'') 

= — — = (5.9) 

Q(l,r) N 

where we call it the filling fraction. According to these fractions, we can introduce the partial 
't Hooft parameters, 

r('''') =g,iV(''^). (5.10) 

Note there are n 't Hooft couplings and filling fractions, but only n — 1 fractions are in- 
dependent since they have to satisfy Yl?=i v^^'^^ = 1 while all the 't Hooft couplings are 
independent. We will show they are related to Coulomb moduli parameters as in the case 
discussed in [651 E] • 

Using these distribution functions, the prepotential and the criticality condition can be 
written as 

= / dxp'^''\x)V{x) 



r=0 ■ 
fc-1 



r=0 
k-l 
i2 ' 



^ P / dxdyp^^^ (x)p(^) (y) log sinh 



x-y 
2 



I dxdypM (x)p(^)(y) log sinh , (5.11) 

r=is \ / 

^V'ix) = Pj dypM(y)coth(^^ 
k-l 



+ E / dyP^'^^y) ^oth ( ^ + ) . (5.12) 



Here P stands for the principal value. To assign the normalization conditions ()5.9p . we then 
introduce the Lagrange multipliers to the original prepotential, 

k~l n , „ ^ 

j-^j^+EEr(''"MT / dxp(")(x) -r('''-M . (5.13) 

This is just the Legendre transformation: P^^'^^ is conjugate to T^''^^. From the variation 
with respect to the density of states, we obtain 



V{x) = 2r y dyp('')(y) log sinh 



2 



fc-i 



+2^E/ '^y/'^'Hy) log sinh + +pa,r)^ ^ g ^a,r)^ 

(5.14) 

which is equivalent to 

v;5:)(x) = r(''^), x^6^^^\ (5.15) 
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In order to solve these conditions, we then introduce the individual resolvents 

N / (r) ' 



(a:) = g,J2coth i^-^ j . (5.16) 

Since all the singularities of the individual ones are found on the real axis, to treat them 
separately, it is convenient to introduce the total resolvent, 

io{x) = J2 ^^"^ - . (5.17) 

Fig. U] shows the singularities of the total resolvent. Their boundary conditions are given by 

— > ibT, u}(x) — > ±kT as x — > ±00. (5.18) 

In the large N limit these resolvents can be also represented with the distribution functions 
as 

oo^^'Hx) = T j dyp^'\y)coih(^^^, (5.19) 

a;(x) = /"(iypW(y)coth('^-^ry (5.20) 

Corresponding to this total resolvent, it is convenient to define the total density of eigenval- 
ues, 

P(x) = |;p(^)(x-^r). (5.21) 

This has supports not only on the real axis but also parallel lines as the total resolvent. Thus 
we can rewrite the resolvent, the prepotential and the saddle point equation in a simple form 
as 

uj{x) = T j dyp{y) coth , (5.22) 

^ = tY^^J dxp^"^ {x)V{x) - T^P j dxdyp{x)p{y) log sinh {^^^ , (5-23) 

^V'ix) = fJ dyp{y) coth + ^r) . (5.24) 

From multi-cut discontinuities of the resolvents we obtain the density of eigenvalues 

/9W(x) = 3— fa;W(x + k)-a;('')(x-ie)') (5.25) 

AniT \ ) 

which has supports on the intervals C'-'''"^. Similarly we have 

^(u;M(x + i6)+^M(x-k)) =P y dypM(y)coth(^^^ 



(5.26) 
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Figure 4: (a) Cuts of the total resolvent on the complex plane. It is cylindrical since there is 
a periodicity in the imaginary direction x ~ x + 27ri. There are totally kn cuts on the plane, 
(b) Expansion around a specific cut. The cycle D^''^) is conjugate to 



2m \ ( 2TTi 
uj \ X -\ — —r + ie \ + uj \x -\ — —r — ie 



Therefore the saddle point equations (|5.24|) can be written as 

(5.27) 

To discuss geometric aspects of the theory, it is convenient to introduce a new function 

y^^){x) = V'{x) -uj{x + = K?'(x). (5.28) 

When we divide the resolvent into regular and singular parts uj{x) = uj^cgix) + casing (x), due 
to the saddle point equation, this function is given by 

yM(x) = -casing (^r+^r^ (5.29) 

We also introduce another function as well as the total resolvents, defined as 

y{x) = kV'ix)-uj{x). (5.30) 

This function has cuts C^^'^^ = + 2r7ri/A;|x G C^^'^'\ I = 1, • • • , n, r = 0, • • • , A; — 1} as shown 
in Fig. m corresponding to all the cuts of the original individual resolvents uj'^'^^ (x). Since the 
resolvent is periodic in the imaginary direction, x ~ x + 27rz, there seems to be the discrete 
Zfc shift symmetry x — > x + 2iii/k. This is quite similar to the discrete Z„ symmetry of SU(n) 
five dimensional theory, proposed in [66j. This similarity suggests the level-rank duality of 
this model. 

Thus the partial 't Hooft coupling is simply given by its contour integral, 

rj^ii^r) = _L I dx Wsing [x + ) = --^ / dx y{x). (5.31) 
47rz Jc(i,r) ^ \ k J 4m Jc(i.r) 

In addition, from (|5.13p we have another contour integral 
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where the contour is the relative cycle represented by the path starting at cxd_|_, and 

going back to oo_ on another sheet after reaching C^^'**) as shown in Fig. HI Indeed this 
relation is analogous to the relation for Seiberg-Witten curve 



ai= f dS, = aj 

lA, dai 



D 



dS. 



(5.33) 



B, 



For our matrix model it is shown in the following discussion that the one-form obtained from 
the spectral curve of the matrix model coincides with Seiberg-Witten differential. 



5.2 Relation to Seiberg-Witten theory 

We now discuss the relation between Seiberg-Witten curve and the matrix model. In the first 
place, the matrix model captures the asymptotic behavior of the combinatorial representation 
of the partition function. The energy functional, which is derived from the asymptotics of 
the partition function [3] , in terms of the profile function 



SaU) = \P dxdyf"{x)f"{y){x - y)^ (log 



(5.34) 



can be rewritten as 
Ea(£>) 



TV 



-P / dxdyp^ - 2 [ dxeix) log TT ^ " 
Ix^y [x-yY J f-J^ 



A 



(5.35) 



up to the perturbative contribution 



{ai - amf log 



Lm 



ai - an 
A 



by identifying 



fix) - - ai\ = q{x). 

1=1 

Then integrating ()5.35p by parts, we have 
Ea(^) = -P / dxdyg' {x)Q'{y)\og{x-y)+2 j dxg'ix)^. 



(5.36) 



(5.37) 



(x -ai) log 



X - ai 
A 



■ (x - ai) 
(5.38) 

This is just the matrix model with the potential (j4.22p with /c = 1 if we identify q'{x) = p{x). 
Therefore analysis of this matrix model is equivalent to that of [3]. But in this section we 
reconsider the result of the gauge theory from the viewpoint of the matrix model (similar 
approach is found in [67j). 

5.2.1 Four dimensional theory 

Let us first consider the most simple case, the four dimensional theory with /c = 1, as a 
preliminary for other generalized theories. This is the case investigated in [3] in detail. We 
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now concentrate on the relation between Seiberg-Witten curve and the spectral curve of the 
matrix model. One can see notations for the four dimensional matrix model in appendix lAl 

To clarify the connection between the matrix model description and the gauge theory 
consequences explicitly, we have to consider the contribution of the dynamical scale, which is 
actually interpreted as the instanton counting parameter. In this case the counting parameter 
should be replaced with A^^""^/)!-^!/'^ because the instanton number becomes fractional |39) . 
In this section we introduce this contribution to the potential functions for both the vector 
multiplet and the (anti) fundamental matter. 

Due to the saddle point equation ()A.7p . we have a regular function on the complex plane, 
except at the infinity, 

Pniz) = A" (e^/^ + = A" (^u; + (5.39) 

This turns out to be a monic polynomial -P„(x) = + • • •, because it is an analytic function 
with the following asymptotic behavior, 

^^n^-uj{z)Y^fi^\ ^ ^n^ ^ ^ CX). (5.40) 

1=1 ^ ^ 

Here w should be the smaller root with the boundary condition as 

A" 

w — > z — ^ oo, (5.41) 

thus we now identify 

w = e-y/^. (5.42) 

Therefore from the hyperelliptic curve ()5.39p we can relate Seiberg-Witten curve to the 
spectral curve of the matrix model, 

dS = z— = logu; dz = - — y(z)dz. (5.43) 

ZTTi w 2m Am 

Note that it is shown in |68[ [57] we have to take the vanishing fraction limit to obtain the 
Coulomb moduli from the matrix model contour integral. This is the essential difference 
between the profile function method and the matrix model description. 

We can deal with the theory with (anti)fundamental matters in a similar way. In this 
case the spectral curve of the matrix model reads 



n Nf 

1=1 f=i 

Thus the regular part is modified as 

PJz) 



(5.45) 
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where we now take into account the contribution from the matter fields, 

Q{z) = Yl{x + mf). (5.46) 
/=i 

As well as the pure gauge theory, we identify w = e~^/^ and Pn{z) is a monic polynomial 
because of its analyticity and asymptotics. 

For k > 1 although we have to consider the multi-matrix model, we can perform almost 
the same approach. In this case we now apply the same resolvent uj and the corresponding 
function y to all the matrices. Due to the slightly modified potential term, we have a regular 
function 

Pn{z) = A" (^e'^y/^ + e~^-^/2^ = A" (^u; + (5.47) 
Again Pn{z) is a monic polynomial Pniz) = + • • •, due to the asymptotic behavior, 

^ ^n^-ku(z) "Q (^^\ . (5.48) 

Identifying the smaller root as 

w = e-*^^/2 — , (5.49) 



z" 

we have a slightly modified relation for the spectral curve 

1 / s , 1 dw 
—y(z)dz = - — z — . 5.50 

We can apply the same spectral curve to each matrix because they are independent of the 
index r of the matrix. Therefore Seiberg-Witten curve for each matrix is given by multiplying 
the original one by 1/k, 

dS yds. (5.51) 

k 

It means Seiberg-Witten curve for the ALE space is decomposed to k curves, and when 
we consider the total contribution from k matrices, we can reproduce the original relation 
for the prepotential and the Coulomb moduli. We can easily see the relation (|5.5ip is also 
obtained from the theory with the matter fields. 

This simple relation between k = 1 and k > 1 theories is also suggested by the result 
in |58j . It is shown that the central charge of the two dimensional conformal field theory, 
corresponding to the theory we discuss here, is given by 

c = fen + ^^^l±^. (5.52) 
k eie2 

In the case of /3 = 1, namely £2 = — ei, since the second part is vanishing, it is natural 
to expect there are k sectors without interacion between them. This is consistent with our 
result discussed above. 
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5.2.2 Five dimensional theory 

We generalize the result for k = 1, which is derived in [66j (see also [3J), to the cases k > 1. 
We now apply a similar method discussed above to the five dimensional case. As well as the 
four dimensional theory, we can have a regular function 

X-^^'/^PkniX) = A"" (^e^^^ + e-y/^') =A'' (^w+ , X = e' . (5.53) 

where we modify the definition y ^ y + kT. In this case it yields a kn^^ monic polynomial 
Pkn{X) = X^"^ + • • •, because of its asymptotics 



1=1 

By setting the boundary condition 



Seiberg-Witten differential for this theory defined on the curve (|5.53|) reads 

dS = -^y{z)dz = -i, log(X) — . (5.56) 

We can see the result of |66] is obtained by setting k = 1. In this case although the curve is 
modified as (|5.53p . the form of the differential itself is not changed. In the four dimensional 
limit i? — > 0, A; sets of n cuts are decoupled because distance between them is Ini/kR as 
shown in Fig. SI A generalization to the theory with matter fields is straightforward. 



6 Further research 

The recent remarkable progress on the four dimensional J\f = 2 theory gives the interesting 
relation to the two dimensional conformal field theory. Thus, according to the results dis- 
cussed above, it is natural to search a two dimensional theoretical counterpart of the orbifold 
theory. In this section we discuss some proposals for such an interesting 2d/4d connection. 

6.1 String theory perspective 

We now discuss a string theory realization of the orbifold theory, and try to find its two 
dimensional description. It is well known that string theory gives us a lot of interesting field 
theoretical consequences, e.g. the holographic approach to a strongly correlated system, the 
method to construct instantons and monopoles and so on. Actually the 2d/4d relation [1] 
can be naturally understood in terms of M-theory. 

Let us consider the string theory realization of A/" = 4 SYM theory on the Taub-NUT 
space [37], which is expected to be related to our case. The Taub-NUT manifold TN^ is 
given by a compactification of the Ai^_i singularity. It approaches to the cylinder x 
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at infinity, and the ALE space is obtained by taking the hmit of i? — )• oo where R stands 
for the radius of S^. We remark although the partition function of the gauge theory on the 
Taub-NUT manifold is not identical to that of the ALE space, we expect this difference does 
not seriously affect the corresponding two dimensional theory because the SU(n) theory on 
both of these manifolds is specified by a flat connection on the boundary as discussed 

later. 

We start with type IIA string theory on TNk x 5^ x with wrapping n D4-branes on 
TA^jfc X . Lifting this to M-theory, we then obtain the compactiflcation of TN^ x x and 
wrapping n M5-branes on TN^ x . By replacing the two-manifold with various Riemann 
surfaces S, we obtain the corresponding M = 2 theories [5]. To study this conflguration we 
go back to type IIA theory with another compactiflcation x x M^. In this case there 
are n D4-branes and k D6-branes wrapping M'^ x and x respectively because the 
circle flbration of the Taub-NUT space has singular points. One can see these D4- and D6- 
branes are intersecting on T^, and thus the chiral fermion is arising from these intersecting 
conflguration. This chiral fermion plays an important role in considering the level-rank 
duality of the system. 

To discuss the two dimensional theory on we then deal with the boundary of the 
four dimensional manifold. By considering the radial quantization near the boundary and a 
wavefunction for the time evolution along x M, we have a Hilbert space with a state 

\p) for each n-dimensional representation 

p: ZkC SU(2) — y SU(n). (6.1) 

Integrating out the flux coming from k D6-branes, we obtain the Chern-Simons term for 

Ics = 2irk [ CS{A). (6.2) 

This means the boundary condition for the D4-brane requires specifying a state of the SU(n) 
Chern-Simons theory at level k living on T^, and thus we have a state for each integrable 
representation of the SU(n)fc WZW model. Note that the same procedure can be performed 
for the SU(/c) theory on k D6-branes, and the diagonal U{1) part is decoupled. Therefore 
we obtain the embedding U(l)nfc x SU(n)fc x SU(A;)„ C \J{nk)i, thus it is easy to see the 
level-rank duality. 

We also remark it is conjectured in [2] that the gauge theory partition function is related 
to the r-function by utilizing the chiral fermion representation. Consequently the generating 
function for instantons on the ALE space of the ADE type would be related to the ADE 
WZW theories on the Seiberg-Witten curve. It is consistent with the intersecting brane 
conflguration because we obtain the Af = 2 theories when we consider the intersecting 
branes on general curves S, which turn out to be Seiberg-Witten/Gaiotto curves [5]. 

We then discuss the connection between the two and four dimensional theory more con- 
cretely. The original proposals of [H] elucidate the explicit relation between the SU(n) 
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partition function and the two dimensional conformal field theory described by Wn alge- 
bra. This relation is extended to the theory with the generalized W algebra |52l [53] [54j 
|55j . which is obtained from the quantum Drinfeld-Sokolov reduction applied to an affine 
Lie algebra SU(n) [56j. It is characterized by a choice of an embedding p : SU(2) SU(n), 
and this embedding p can be labeled by a partition of n, or equivalently a Young dia- 
gram Y. For example, it reproduces well-known algebras, VF(SU(n), [1, • • • , 1]) = SU(n), 
W^(SU(n), [N]) = Wn, etc. Such an embedding has been discussed in the orbifold case ()6.ip . 
Indeed we consider the embedding of the finite subgroup of SU(2) for the orbifold theory, 
and the decomposition into the irreducible representations of corresponds to that of the 
gauge group SU(n) — > SU(no) x • • • x SU(nfc_i) with no + • • • + rif^^i = n. This shows 
that this embedding is also labeled by a partition of n. Therefore it is expected that the 
two dimensional conformal field theory, corresponding to the SU(n) gauge theory on the 
ALE space C^/Z^, is related to the generalized W algebra W{S\]{n), [no, • • • ,nfc_i]). This 
connection is still speculation, and should be investigated in detail. 

6.2 g-deformed CFT 

It is shown in that g-deformed CFT, which is described by g-Virasoro [70] and q-W 
algebra [7l], corresponds to the g-deformed five dimensional partition function of the gauge 
theory, and its matrix model description is also proposed in [72] by clarifying a connection 
with g-Virasoro algebra. This matrix model, which we now call g-Virasoro matrix model, is 
apparently different from the trigonometric one |26[ [271 [28] , but it is worth considering the 
root of unity limit of such a Q'-Virasoro related models. The root of unity limit of g-Virasoro 
algebra is investigated in [73] while that of q-W algebra is not yet well known. 

In this subsection let us focus on the matrix measure part of the q-Virasoro matrix model 
[72] . It is just given by the two-parameter deformed Vandermonde determinant, which is 
closely related to the Macdonald polynomial [74] , 

We can obtain the usual Vandermonde determinant by taking the limit g — )■ 1 with t = q^, 

^ n (i - f")' - iK-^ - (6-4) 

This corresponds to the limit to get the Jack polynomial from the Macdonald polynomial [73] . 
On the other hand, performing the Uglov condition, q — )• ojq, t = uq^ with (3 = + 1 = 1 
(mod k), then taking q — > 1, it becomes 

7 



^ \ ""jj i<j 
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i<j 

(6.5) 

The matrix measure (|6.5p is expected to be related to the corresponding one (|4.14|) . Indeed 
they are equivalent for k = 1. Although this is written in terms of only one matrix while 
the model discussed in the previous section is not, we can see the Vandermonde part in (j6.5p 
is also found in ()4.14p . We just expect the interaction part in ()4.14p is encoded into the 
non-singular part of ()6.5p . 

Such a difference seems to be related to that found in the usual matrix model for the 
four dimensional theory, the Penner-type [HI [161 El ttS] , interpreted as the integral rep- 
resentation of the Liouville correlator and the other one just given by the asymptotics of 
the combinatorial expression of the partition function [26t [271 128] , Indeed they are differ- 
ent because the conformal symmetry is manifest for the former one, but the latter one. 
However it is not enough yet to understand the difference between the Macdonald-type and 
trigonometric-type deformed Vandermonde determinant. This difference is expected to be 
understood, for example, in terms of string theory. 

We now remark the polynomial, called the Uglov polynomial in |42] . is given by taking 
this limit from the Macdonald polynomial in order to describe the spin Calogero-Sutherland 
model [H]. 

7 Summary and discussion 

In this paper we have performed some extensions of the partition function and its matrix 
model description for the four dimensional J\f = 2 gauge theory to the orbifold theory. 
We have shown that the orbifold projection, extracting the F-invariant sector of the Young 
diagram, is automatically performed by taking the root of unity limit q exp{2TTi/k) of 
the g-deformed partition function while we have to take g — > 1 to get the four dimensional 
theory on M^. For such an orbifold partition function of SU(n) theory with F = Z^, it is 
convenient to divide n-tuple to A;n-tuple partitions. As a result, we have obtained the multi- 
matrix model by considering the asymptotic behavior of the combinatorial representation. 
This matrix model at the large limit has been analysed in detail, which is equivalent 
to studying the limit shape of the Young diagram, and then we have seen Seiberg-Witten 
curve is obtained as the spectral curve of the matrix model. We have also discussed the 
corresponding two dimensional theory of the four dimensional orbifold theory. Focusing on 
the embedding p : C SU(2) — )• SU(n), which characterizes the decomposition of SU(n) 
gauge group to the irreducible representations of Z^, we have suggested the generalized W 
algebra appears in the two dimensional theory. 

We now comment on some possibilities of extension beyond this study. We hope our 
study becomes a step for understanding of M-theory itself. Actually the emergence of S\J{n)k 
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symmetry in the two dimensional theory as a counterpart of the orbifold theory C'^ /'Lj. is 
analogous to the ABJM theory because the level of Chern-Simons theory in the ABJM the- 
ory is directly related to the background manifold C^/Z^ on which M2-branes are located. 
Indeed, in both cases, the level k is interpreted as the degree of singularity of the complemen- 
tary manifold of the corresponding world- volume theory. It is interesting to study a relation 
between M2- and M5-branes from this point of view. We are also interested in the level-rank 
duality of our model. Such a duality can be found in the two dimensional SU(n)fc, or SU(A;)„ 
WZW theory, which is closely related to the four dimensional orbifold theory. It is expected 
that this duality plays an important role in understanding some aspects of M-theory. 

It is natural to consider some applications to related topics. One of them is the three 
dimensional duality, which is a recent hot topic on this subject [73 EUJ [771 ESI ESI [Sni EI]- 
The g-parameter plays a similar role in such a theory, so that it is interesting to study the 
singular limit of the g-parameter, i.e. the root of unity limit. Second is a relation to the 
quantized integrable models. It is shown in [82j that when we consider generic 0-parameters, 
ei + €2 7^ 0, correction to the prepotential can be interpreted as a quantization effect of the 
corresponding integrable model. Searching an integrable model, which corresponds to the 
orbifold theory, would be also interesting. 
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A Four dimensional matrix model 

In this appendix we summarize the results of the four dimensional limit of the matrix model 
in order to fix our notations. Let us mainly consider the fc = 1 theory, and shortly comment 
on simple generalizations to the cases of A; > 1. 
The partition function for this case is given by 



and the saddle point equation is obtained from its differentiation with respect to each eigen- 
value, 




(A.l) 




(A.2) 
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This is also given by the extremal condition of the effective potential defined as 

N 

Ves{xi) = V{xi) - 2gs ^ log(xi - Xj). (A.3) 

We then introduce the resolvent for this model. By taking the large N limit, it can be 
given by the integral representation, 

u(x) = T [ dy-^. (A.4) 
J x-y 

Its asymptotic behavior yields 

uj(x) — >■ — , X — > oo. (A. 5) 

Here the density of states is obtained from the discontinuities of the resolvent, 

Thus the saddle point equation can be also written in the following form, which is convenient 
to discuss its analytic property, 

V'{x)=uj{x + ie) + uj{x-ie). (A.7) 

On the other hand, wc have another convenient form to treat the saddle point equation, 
which is called the loop equation, given by 

- V'{xf + R{x) = (A.8) 

where we denote 

y(x) = V'{x) - 2uj{x) = -2uJsuigix), 
i=l 

It is obtained from the saddle point equation by multiplying l/(x — Xi) and taking their 
summation and the large N limit. This is interpreted as the hyperelliptic curve which is 
given by resolving the singular form, 

y'^ix) - V'{xf = 0. (A.IO) 

The genus of the Riemann surface is directly related to the number of cuts of the corre- 
sponding resolvent. The filling fraction, or the partial 't Hooft coupling, is simply given by 
the contour integral on the hyperelliptic curve 

7] = ^ ^ dx LUsingix) = ^ dx y(x). (A.ll) 
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For A; > 1 we can perform the same thing as well as the case of A; = 1 since there is 

no interaction between k matrices in the four dimensional limit with (3 = 1. The /c-matrix 
partition function is given by 

k-l 

Z=Y[Z^'^\ (A.12) 

r=0 

= /n^n(xr -xf )^e-^S-^(^^^'). (A.13) 
i=i i<:j 

We can provide all the matrices with the same density of states, resolvents and hyperelliptic 
curves etc, e.g. y^'^\z), and the total contribution to the partition function is given by the 
simple summation of them. On the other hand, when we consider the case /3 ^ 1, we have to 
take into account interaction between matrices. It is a quite interesting situation, and will 
be investigated in a future work. 
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